Correlated electrons in Fe-As compounds: a quantum chemical perspective 
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State-of-the-art quantum chemical methods are applied to the study of the multiorbital correlated 
electronic structure of a Fe-As compound, the recently discovered LiFeAs. Our calculations predict 
a high-spin, S~2, ground-state configuration for the Fe ions, which shows that the on-site Coulomb 
interactions are substantial. Also, orbital degeneracy in the (xz,yz) sector and a three-quarter 
filling of these levels suggest the presence of strong fluctuations and are compatible with a low 
metallic conductivity in the normal state. The lowest electron-removal states have As 4p character, 
in analogy to the ligand hole states in p-type cuprate superconductors. 



Electron correlations in transition-metal (TM) solid- 
state compounds give rise to a variety of less conven- 
tional phenomena. Superconductivity in layered copper 
oxides at temperatures as high as 100 K, for example, and 
the pairing mechanism in these systems are believed to 
involve strong correlation effects among the valence Cu 
3d electrons. The discovery of the Fe 2+ -As 3_ supercon- 
ducting compounds [l[ is the latest surprise in the field of 
(i-electron systems. Also in this case, the experimental 
data indicate unconventional superconductivity. "Bad- 
metal" conductivity in the normal state, a small carrier 
density, a relatively small in-plane coherence length, and 
Uemura scaling in the muon spin relaxation spectra 
are all reminiscent of cuprate superconductors. 

How strong correlations are in iron pnictides is an im- 
portant issue. In a first approximation, the charge distri- 
bution within the Fe 3d levels and the spin state depend 
on the intra-orbital Coulomb interaction, the so-called 
Hubbard U. A U value much larger than the crystal-field 
splittings will favor a high-spin arrangement of the six Fe 
3c? electrons, while a low-spin, closed-shell configuration 
is expected for small values of U. The picture complicates 
when inter-orbital Coulomb and exchange interactions 
are considered. Additionally, correlation effects related 
to ligand p to TM d charge transfer excitations may be 
important too, as discussed for the case of another 3d 6 
system, the cobalt oxide perovskite LaCo03 

We investigate here the electronic structure of a Fe-As 
compound, the recently discovered LiFeAs Ab ini- 
tio, wave-function-based methods from modern quantum 
chemistry are used in our study. We characterize the 
ground-state electron configuration for the undoped case 
and provide new insight into the nature of doped car- 
riers. Our results lend credence to the view [5|, |6j that 
correlations are moderate to strong in Fe pnictides. 

LiFeAs has a tetragonal crystal sructure, with the 
VA/nmm space group Different from other Fe- 

As compounds, it exhibits superconductivity at ambi- 
ent pressure without chemical doping, with T c k, 18 K. 
The common feature of the Fe pnictide superconduc- 
tors is the Fe2As2 network of FeAs4 tetrahedra. Near- 
est neighbor (NN) FeAs4 units share edges, while next- 
nearest-neighbor (NNN) tetrahedra share their corners. 



In LiFeAs, the Fe2As2 layers are separated from each 
other by double layers of Li ions. 

The first step in our study is a ground-state restricted- 
Hartree-Fock (RHF) calculation for the periodic crystal. 
This calculation was carried out with the CRYSTAL pack- 
age hf. We employed the lattice parameters reported in 
Ref. [4[ and Gaussian-type, all-electron basis sets. Basis 
sets of double-zeta quality from Towler's CRYSTAL data 
basis were applied for the As and Li ions [1, 0]. For 
Fe, we used a basis set of triple-zeta quality, with s and p 
functions from Towler's data basis [8J] and the d functions 
developed by Seijo et al. [loj ]. 

The periodic RHF calculation yields a finite gap at the 
Fermi level. For the RHF wave function, the xz and yz 
components are the highest among the Fe 3d levels and 
unoccupied; the other Fe 3d orbitals are doubly occupied. 
We employ a reference system having the x and y axes 
rotated by 45° with respect to the a and b coordinates 
of the P4/nmra space group, such that the As NN's of a 
given Fe site are situated either in the xz or yz plane. 

On-site and inter-site correlation effects are in- 
vestigated in direct space by means of multi- 
configuration complete-active-space self-consistent-field 
(CASSCF) and multireference configuration-interaction 
(MRCI) calculations. The CASSCF wave function is con- 
structed as a full configuration-interaction (CI) expan- 
sion within a limited set of "active" orbitals [llj, i.e., all 
possible occupations are allowed for these active orbitals. 
In the present study, the active orbital set contains all 3d 
orbitals at a given number of Fe sites. Not only the CI 
coefficients but also the orbitals are variationally opti- 
mized in CASSCF, which makes this method quite flex- 
ible. MRCI wave functions are further constructed by 
adding single and double excitations from the Fe 3s, 3p, 
3d and As 4s, Ap orbitals on top of the reference CASSCF 
wave function, which is referred to as SD-MRCI 
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The quantum chemical computations are performed on 
a finite cluster C including nine FeAs4 tetrahedra - a 
"central" FeAs 4 unit plus four NN and four NNN tetrahe- 
dra - and 16 Li neighbors of the As ions of the "central" 
unit. The orbital basis entering the correlation treat- 
ment is a set of projected RHF Wannier functions: lo- 
calized Wannier orbitals (WO's) are first obtained with 
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the Wannier-Boys localization module [l2j of the CRYS- 
TAL package and subsequently projected onto the set of 
Gaussian basis functions associated with the atomic sites 
of C Projected As Ap x and Fe 3d xz WO's, for ex- 

ample, are plotted in Fig. 1. Moreover, the RHF data 
is used to generate an embedding potential for the nine- 
tetrahedra fragment C. This potential is obtained from 
the Fock operator in the RHF calculation [l3( and models 
the surroundings of the finite cluster, i.e., the remaining 
of the crystalline lattice. It is added to the one-electron 
Hamiltonian in the subsequent CASSCF/MRCI calcula- 
tions via an interface program developed in our labora- 



tory [1J|. The CASSCF and MRCI investigations are 



carried out with the MOLPRO program [if 

In a first set of CASSCF calculations, a number of 
nine sites are explicitly correlated, those of the reference 
FeAs4 tetrahedron and the four NN Fe ions. This group 
of nine sites form the active region of the cluster, which 
we denote as Ca- The other ions in C, i.e., four Fe NNN's, 
12 As, and 16 Li sites, form a buffer region Cb, whose 
role is to ensure an accurate representation of the tails of 
the WO's centered in the active part Ca- For our choice 
of Cb, the norms of the projected WO's at the central 
tetrahedron and NN plus NNN Fe sites are not lower than 
99.5% of the original crystal WO's. While the occupied 
orbitals in the buffer zone are frozen, orbitals centered 
at sites in the Ca region (and their tails in Cb) are al- 
lowed to relax and polarize in the CASSCF study. For 
this first set of CASSCF calculations, the active orbital 
space consists of 25 Fe 3c? orbitals, five at each active Fe 
site. Due to the large size of our CAS, i.e., 30 electrons 
in 25 orbitals, we restrict our calculations to the case 
of high-spin (i.e., ferromagnetic) couplings among neigh- 
boring Fe ions, although the experiments indicate antifer- 
romagnetic inter-site interactions. It is known, however, 
that the local charge distribution at a given TM site does 
not depend on the nature of the inter-site d-d magnetic 
couplings. We observed this in the case of the Cu ox- 
ide superconducting compounds, for example. Also, a 
closed-shell representation of the Fe ions in the rest of 
the crystal is acceptable, given the fact that the Fe NN's 
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FIG. 1: a) As Ap x and b) Fe 3d xz RHF WO's after projection 
onto the finite cluster, see text. The tails at nearby sites are 
very small. For the Fe 3d xz and 3d yz WO's, the lobes pointing 
towards the NN As sites are less extended. 



are treated by CASSCF, on the same footing with the 
reference Fe site. 

The CASSCF calculations show that at each Fe site the 
3d electrons are coupled into quintet states. In contrast 
to the RHF results [16| , the xz and yz components are the 
lowest in energy, such that the sixth electron is accomo- 
dated into these levels. The ground-state wave function 
is thus doubly degenerate, d xzyz d xv d\ z2 _ r2 d 1 x2 _ y2 . The 
first excited state is also a quintet, at 0.25 eV higher en- 
ergy, and corresponds to a {xz,yz)—+xy transition [171 ] . 

The fact that the relative energies of the Fe 3d lev- 
els are not consistent with simple considerations based 
on ligand-field theory for tetrahedral coordination, which 
predicts that the e g levels are lower than the ti g states, 
was pointed out before, see, e.g., Refs. @,[I3]. This seems 
to be related to the distortion of the As tetrahedra, which 
are squeezed in the z direction, and the presence of direct, 
nearest- neigbor Fe d-d orbital overlap 0, [l8| . 

Our finding of a high-spin (HS), 5 — 2 ground-state 
configuration agrees with the results of dynamical mean- 
field theory (DMFT 19]) investigations in the moderate 
to strong coupling regime by Haule et al. [j| and Craco et 
al. [6J]. Also, a three-quarter filling of the degenerate d xz 
and d yz bands as found in our calculations is compatible 
with the low metallic conductivity in the normal state of 
these systems. Systems where a pair of degenerate or- 
bitals accomodates one electron or one hole often display 
very rich physics, involving couplings among the charge, 
lattice, and spin degrees of freedom. The structural tran- 
sition at about 150 K in some of the Fe-As compounds 
"2(il 21 1 might occur such that the degeneracy of the xz 



and yz orbitals is lifted. This issue remains to be inves- 
tigated in future work. 

It would be instructive to determine the relative ener- 
gies of states involving low-spin couplings at a given Fe 
site. However, for a CAS with five Fe ions and 25 orbitals, 
such investigations are quite difficult. Since the on-site 
interactions are much larger than the inter-site d-d spin 
couplings, the states related to low-spin configurations 
at a given site are among the highest in a multitude of 
low-spin excited states. Identifying and optimizing those 
states is a very tedious task. In order to access those 
states, we reduce then the orbital space in the CASSCF 
calculations to the set of five 3d functions at the central 
Fe site. For each of the Fe NN's, the 3d electrons are 
"forced" into a t\ g closed-shell configuration. 

With this choice of the CAS, the lowest intermediate- 
spin (IS), 5 = 1 and low-spin (LS), 5 = states require 
excitation energies of 1.91 and 2.34 eV, respectively, with 
respect to the HS ground-state (see Table I). The lowest 
5 = state corresponds to a t® g orbital occupation and 
the lowest 5 = 1 state corresponds to a t\ g e g configura- 
tion, more precisely d xy d\ zyz d\ z2 _ r2 . 

On top of the CASSCF wave functions, we further per- 
formed MRCI calculations with single and double exci- 
tations from the Fe 3s, 3p, 3d and As 4s, Ap orbitals at 
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TABLE I: Relative energies for the HS, IS, and LS configura- 
tions by CASSCF and SD-MRCI calculations, see text. The 
HS state is always the lowest. The energy of the HS state in 
the SD-MRCI calculation was taken as reference. 



Relative energy (eV) 


HS 


IS 


LS 


CASSCF 


13.16 


15.07 


15.50 


SD-MRCI; Fe 3s,3p,3d, As 4s,4p 





1.30 


1.74 



the central tetrahedron. The SD-MRCI treatment de- 
creases the HS-IS and HS-LS energy splittings to 1.30 
and 1.74 eV, respectively (lowest line in Table I). Now, 
the IS-LS energy difference, for example, can be used 
to extract information on the magnitude of parameters 
such as the intra-orbital and inter-orbital Coulomb in- 
teractions U and U' and the inter-orbital exchange cou- 
pling Jh. The t\ g e g — t% g energy difference can be ex- 
pressed as AE = U + Jh — U' — Acf- For the sake of 
simplicity, we assume that the Coulomb repulsion terms 



between different t 2g and e g orbitals have all the same 
value. A unique value is also assumed for the t 2g —e g inter- 
orbital exchange couplings [13] • To determine the crystal- 
field splitting Acf between the (xz,yz) and 3z 2 — r 2 
components, we perform SD-MRCI calculations for the 
(xz,yz) — > 3z 2 — r 2 excitation energy in the HS, S — 2 con- 
figuration. By SD-MRCI, this splitting is 0.71 eV. Since 
AE = 0.44 eV, see the lowest line in Table I, it follows 
that U + J H -U' = 1.15. Using the relation U' = U-2J H , 
we find [/-[/'« 0.8 eV and J H w0.4 eV. 

As expected, U — U' does not depend on the electron 

configuration at the NN Fe sites: if in the CASSCF 

and SD-MRCI calculations for the central tetrahedron 

we adopt a eld 2 .. , configuration at the NN Fe sites, as 
y x y ^ I — I 

found in the periodic RHF calculation [16| , the difference 
between U and U' remains 0.8 eV. Our results provide a 
lower limit for the value of the Hubbard U. We note that 
estimates based on density-functional (DF) calculations 
strongly depend on the type and size of the Wannier-like 
orbital basis. Constrained DF computations by Anisimov 
et al. [H using a WO basis restricted to the Fe 3d or- 
bitals yield [7 = 0.55 and Jjy = 0.5 eV. With an extended 
orbital basis including As Ap functions, it was found that 
U = 3-j-4 and = 0.8 [23j. Constrained random-phase- 
approximation (RPA) calculations 24] on top of the DF 
data lead to U = 2.2 = 3.3 and J H = 0.3 = 0.6. Values of 4 
eV were used for U in DMFT investigations by Haule et 
al. [H and Craco et al. A value of [7 = 0.3 eV was 
employed by Korshunov and Eremin [25] ] in RPA calcu- 
lations for the spin response in the normal state of Fe 
pnictide compounds. 

In Ref . [2^1 , Kroll et al. found that model-Hamiltonian 
multiplet calculations with [7=1.5 and Jh = 0.8 = 0.9 re- 
produce well the x-ray absorption Fe T^^-edge spectra. 
Regarding the crystal-field splittings, these authors used 



a value of 0.25 eV for the t 2g —e g energy separation, sub- 
stantially lower than our SD-MRCI result of 0.71 eV for 
the (xz,yz) — > 3z 2 — r 2 excitation energy. In order to re- 
produce the experimental spectra, a larger value for the 
crystal-field splitting would imply a larger value for U in 
the model-Hamiltonian calculations. 

That the Coulomb interactions are substantial in these 
systems is best illustrated by the nature of the lowest 
electron-removal states. Our calculations show that for 
these states the additional holes populate the ligand p 
levels, which resembles the situation in layered cuprates. 

In a simple picture, the formation of oxygen 2p hole 
states in p-type cuprates is due to the large Coulomb re- 
pulsion at the Cu 3d 9 sites [233 : in order to minimize 
the interaction with the Cu 3d holes, extra holes in the 
doped system enter the O 2p levels. Ab intio quantum 
chemical calculations show indeed that the first (N — l) 
states have O 2p character in copper oxides, see, for ex- 
ample, the analysis in Ref. [28]. The situation is quite 
similar in LiFeAs, where we find that for the lowest ion- 
ized state the extra hole is accommodated into the As 
Ap orbitals. Some details are, however, different. Due to 
the xz and yz degeneracy, the ligand hole is distributed 
with equal probability over two linear combinations of As 
Ap orbitals, i.e., a "+" combination of slightly tilted p x 
functions at As sites in the xz plane, see Fig. 2, and a 
combination of tilted p y functions at As sites in the yz 
plane. Not only the on-site Coulomb repulsion, but also 
the inter-site correlations are effective: in the CASSCF 
wave function \^>) for the lowest ionized state, the largest 
weight is aquired by configurations where the two holes in 
the (xz,yz) sector reside in pairs of "orthogonal" orbitals, 
either d yz and p x i2=Pxi+Px2 or d xz and p v 34 =p y3 +p y i ■ 
Those are the first two terms in the CASSCF expan- 
sion I*} = 0.63|j£ 12 P 2 y3i d 2 xz dl z } + 0.63^1^1^1^) + 
0.33\p 2 xl2 pl 34 dl z dl z }-0M\p 1 xl2 pl 3i d 2 xz dl z )+... 1 where the 
subscripts 1,2 and 3,4 refer to As sites in the xz and yz 
planes, respectively, and the HS-coupled electrons in the 
Fe xy, x 2 —y 2 , and z 2 orbitals are omitted. Also dif- 
ferent from cuprates, the ligand p and TM d holes are 
HS-coupled for the lowest (N—l) state, with a total spin 
5 = 5/2. Configurations where both (xz,yz)-\ike holes re- 
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FIG. 2: a) Linear combination of Ap hole orbitals at As sites 
in the xz plane for the lowest electron-removal quasiparticle 
state, see text, b) Fe 3d xz orbital for the lowest (N—l) state. 
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side at the TM site contribute as well to the (N— 1) wave 
function, see the third term in the above expansion, as 
also found in cuprates for the two holes of x 2 — y 2 sym- 
metry (2R [2^|. In the Fe-As system, however, the two 
holes do not occupy the same d orbital. The fourth term 
in the expression of \^f) refers to a configuration where 
both holes have As 4p character. 

The CASSCF results for the lowest (iV-1) state |*) 
were obtained with a CAS that contains five Fe d and two 
As p orbitals. The 3d orbitals of the Fe NN's were kept 
again in a t\ g configuration. Also of interest is the na- 
ture of the next higher-lying ionized state. That involves 
a hole in a composite As 4p orbital of x 2 — y 2 symme- 
try extending over all four ligands around the TM site, 
like the Zhang- Rice hole in cuprates [27]], LS-coupled to 
the hole in the Fe d x 2_ y 2 orbital. The effective occupa- 
tion numbers at the Fe site are d xzyz d xy d\ z2 _ r 2d 1 x2 _ y2 , 
with a HS on-site coupling, such that the total spin for 
the FeAs 4 tetrahedron is 5 = 3/2. By CASSCF, the en- 
ergy separation between the two hole states is 0.59 eV, 
where the active space for the higher state Y$>') contains 
six orbitals. MRCI calculations with single and double 
excitations from all Fe 3c? and As Ap orbitals at the FeAs4 
tetrahedron yield a splitting of 0.75 eV. We also inves- 
tigated the nature of the lowest electron-addition states. 
Our calculations show they have Fe (xz,yz) character. 

Ab initio quantum chemical calculations for determin- 
ing the dispersion of the (A=f1) quasiparticle bands and 
the Fermi-surface topology of Fe pnictide systems are 
left for future work. A delicate issue is in this context 
the treatment of renormalization effects due to inter-site 
spin interactions and spin-polaron physics. In cuprates, 
we found 28j that such effects lead to a renormalization 
of the NN hoppings by a factor of 4. In Fe pnictides, 
the experiments indicate effective mass renormalization 
factors of about 2 as compared to the DF data Hi EI. 

To summarize, we apply multiconfiguration CASSCF 
and multireference CI methods to the study of the mul- 
tiorbital correlated electronic structure of a Fe-As com- 
pound, the recently discovered LiFeAs. On-site and Fe- 
As inter-site correlation effects are treated on equal foot- 
ing in our approach. Our calculations predict a HS 
ground-state configuration for the Fe ions, in agreement 
with DMFT calculations Q for systems from the same 
family and simulations of the x-ray absorption spectra 
[201 . The lowest electron- removal quasiparticle states 
have As 4p character, in analogy to the ligand hole states 
in p-type high-T cuprate superconductors. The results 
indicate that the on-site Coulomb interactions are sub- 
stantial. We find that U — U' « 0.8 eV, which provides a 
lower bound for U, and the Hund coupling constant Jh 
is about 0.4 eV. Also, orbital degeneracy in the (xz,yz) 
sector and a three-quarter filling of these levels suggest 
the presence of strong fluctuations and are compatible 
with a "bad-metal" conductivity in the normal state. 
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